Unavoidable strong magnetic fields in the early Universe 
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In this letter we show that the Universe is already strongly magnetised at very early epochs during 
cosmic evolution. Our calculations are based on the efficient amplification of weak magnetic seed 
fields, which are unavoidably present in the early Universe, by the turbulent small-scale dynamo 
(SSD). We identify at least one epoch during the radiation dominated regime where all the necessities 
for the SSD to work are fulfilled. Hence, at scales of l c ~ 0.3 pc, the comoving field strength due to 
this mechanism will be -Bo ~ 0.35s 1 ' 2 nG at the present time, where e is the saturation efficiency. 



INTRODUCTION 

Magnetic fields of strengths of order a few pG have 
been observed in high redshift as well as present day 
galaxies, and in galaxy clusters and superclusters (see 
e.g. refs. [l|-|j] respectively). There is also evidence for 
strong extragalactic magnetic fields coming from 7-ray 
observations. These studies place a lower bound on in- 
tergalactic magnetic fields at 3 x 10~ 16 G [5], although 
plasma effects may complicate the propagation of electro- 
magnetic cascades [a, 0] ■ Theoretically, magnetic fields 
are very likely to have been generated at some level in the 
early Universe through a variety of mechanisms. On the 
largest scales, magnetic fields can be generated during 
inflation giving today Bq ~ 10 -25 — 10 _1 nG on a scale 
of 1 Mpc [8[. However such mechanisms require some 
modification to the Maxwell theory in order to break 
its conformal invariance. Without such modifications, 
fields of present strengths B$ ~ 1CU 20 and ~ 1CU 11 nG 
could have been generated at the electroweak and QCD 
phase transitions respectively [9(]. At later times, mag- 
netic fields could have been generated through the gen- 
eration of vorticity in the plasma of the early Universe 
10H13J (originally pro posed by Harrison [14|). The mech- 
anism in refs. [id . Ill| is very natural, since vorticity in 
the plasma is unavoidably generated in the late radia- 
tion era through the non-linear couplings of first order 
density perturbations. The seed field generated here is 
Bq ~ 10~ 20 nG. Therefore, in order to explain observa- 
tions, some amplification of the seed fields must have 
occurred at some point in the history of the Universe. 
A popular mechanism for such amplification is known as 
the dynamo mechanism. 

The dynamo mechanism has two classes (see ref. [15[ 
for a review). The large-scale dynamo converts kinetic 
energy on large scales into magnetic energy. This mech- 
anism can act if the conducting fluid flow is highly he- 
lical, inhomogeneous or anisotropic, the typical example 
being the differential rotation of galaxies. Whereas the 
SSD, which operates on much smaller scales, converts 



kinetic energy from turbulent motion into magnetic en- 
ergy. Magnetic field lines which are frozen in the con- 
ducting plasma will be stretched by the random motion 
of fluid elements leading to field amplification. The SSD 
can take effect for nonhelical, homogeneous and isotropic 
turbulent flow. This mechanism has so far only been im- 
plemented to the formation of the first stars and galaxies 
in the matter dominated Universe [laj (see ref. 17 1 for an 
early discussion on this subject). In this paper we inves- 
tigate the SSD amplification of a seed field in the radi- 
ation era. We demonstrate that the conditions are right 
for efficient dynamo amplification leading to large mag- 
netic fields which could explain observations and have an 
impact on early structure formation. 



TURBULENT FLUID VELOCITY IN THE 
RADIATION ERA 

To explain the formation of the Large Scale Structure 
(LSS) in the Universe observed today, a primordial den- 
sity perturbation of magnitude Sp/p ~ 1CP 5 is required 
at the time of matter-radiation equality. The primordial 
density perturbation is thought to have been generated at 
a much earlier time and therefore must be present in the 
very early Universe during the radiation dominated (RD) 
era. Cosmic inflation provides the most compelling expla- 
nation for the origin of the primordial curvature pertur- 
bation. A period of quasi-de Sitter inflation generically 
predicts a spatially flat Universe overlaid with a primor- 
dial spectrum of adiabatic, Gaussian and (nearly) scale- 
invariant super horizon perturbations. Inflation generates 
scalar (density) and tensor (gravitational waves) modes. 
Vector modes, at linear level, are not sourced by inflation. 
The theory is in excellent agreement with observations of 
the LSS and CMB anisotropies 
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Well before horizon entry, the primordial curvature 
perturbation, which determines the gravitational poten- 
tial $, remains constant and given by the initial condition 
$0. In the Newtonian gauge, the equation of motion for 



the gravitational potential is given by [20j 



(1) 



where p = wp, % = a' /a and a prime denotes differenti- 
ation with respect to conformal time r/. In the RD era 
w = 1/3, the solution to the above for the Fourier modes 
reads 



$( fe)?7 ) = ^M$ (fc), 



(2) 



where y = krj/^/3 and ji(y) = s'my/y 2 — cosy/y is the 
first spherical Bessel function. As a Fourier mode of the 
gravitational potential re-enters the horizon during RD it 
begins to oscillate with an amplitude decreasing as 1/y 2 . 
The initial conditions, which are probed by observa- 
tions, are given by the two-point correlation function 



($o(fti)$S(ft 2 )) = (2nfP <s> (k)S 3 (k 1 - fc 2 ), (3) 



where 
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From the Planck results A^(fc ) = 2.215 x 10~ 9 and 
n s ~ 0.96 for the pivot scale fc = 0.05 Mpc -1 [19]. Mod- 
els of inflation are constructed in order to generate this 
amplitude and slight scale-dependence. 

Perturbations in the fluid 3- velocity field are generated 
by the density perturbations. At first order in density 
perturbations, the fluid velocity perturbation is purely ir- 
rotational (curl-free) i.e. compressive turbulence, whose 
Fourier modes arc 1201 



Vi(k,T)) 



2H 2 



[&(k,TJ)+H$(k,ri)] 



With the solution in eq. ([2]) we find 
i3v3fc; 



Vi(k,Tj) = -- 



siny-2ji(y)]$ (fc). 



(5) 
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The perturbation in the velocity field oscillates with the 
density perturbation, however it has a term which does 
not decay with the expansion. The spectrum of velocity 
perturbations can be defined by the two-point correlation 
function 



(v i (k 1 ,ri)v*(k2,v)) = (2^) 3 P„(fc),5 3 (fe 1 
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where V v (k) = (k 3 /2-K 2 )P v (k). Considering eq. | 
find 

27 
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therefore 
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For subhorizon scales y, k/H ^> 1 the spectrum oscil- 
lates rapidly, so we can average the spectrum over many 
oscillations. For a scale-invariant primordial spectrum 
n s = 1 we find 
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V v {k) ~ ^A 2 j(fc ). 



(10) 



Thus, on subhorizon scales the spectrum of velocity per- 
turbations becomes scale-invariant. The typical value of 
the velocity perturbation is therefore 

v ~ ^Vv(k) ~ y/A^(ko) ~ 5 x lO" 5 , (11) 

where natural units are used such that c = 1. Hence, the 
turbulent velocity spectrum generated by first order den- 
sity perturbations is isotropic, homogeneous, Gaussian 
and to a good approximation scale- invariant v(l p ) ex it 
with •& = on a physical scale l p . 
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FIG. 1. The spectrum of velocity perturbations at temper- 
ature T = T f ~ 19.6 keV. The comoving length l c ~ 0.27 pc 
(black-dashed) is well within the oscillatory regime. (Color 
online) 



THE SMALL-SCALE DYNAMO 

In the SSD mechanism, kinetic energy from turbulent 
motion is converted to magnetic energy. There is a large 
number of numerical works in the literature which have 
demonstrated this effect unambiguously for a number of 
settings 2ll424j. Generically, the fluctuating component 
of the magnetic field grows exponentially B oc exp(ri), 
where the growth rate T depends on the kinetic Reynolds 
number R e . Analytically, the Kazantsev theory (follow- 
ing the formalism by refs. 15s, 25]) was developed in order 
to study the evolution of magnetic fields in a conducting 
plasma containing turbulent motion 26]. For random 
motions correlated on a scale l p with velocity v <xlt, un- 
der certain conditions it can be shown that the growth 



rate is given by 27 | 



(163-304tf)v 
60L 



R (i-#)/(i+4) _ 



(12) 



This result is obtained in the large Prandtl number 
P m = R m /R e limit, where the Prandtl number is a mea- 
sure of the relative importance of the magnetic and ki- 
netic diffusion. Here, the R m is the magnetic Reynolds 
number. We note that the above analytical result was 
obtained by modelling the turbulent velocity spectrum 
so that it behaves as Kolmogorov or Burgers turbulence 
for scaling index $ = 1/3, 1/2 respectively [27j (see con- 
clusions for further discussion). 

There is a critical magnetic Reynolds number for which 
R m > R™ and T > 0, and the SSD takes effect. The 
magnetic Reynolds number is given by R m — Attlpva 
where a is the plasma conductivity. The local kinetic 
Reynolds number R e characterises the relative impor- 
tance of the fluid advective and dissipative terms in the 
Euler equation. On some physical length scale l p , it is 
given by 



R e {l p ) 



v(l p )l p 



v(lp)l c 
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(13) 



where l c is a comoving length. The shear viscos- 
ity T) s = {9v,i/^'9*)lp 1 IS determined by the particles of 
longest mean free path {mfp) lp' 1 , which are either neu- 
trinos or photons in the radiation dominated era. Here, 
5* and g vn are the total and component number of ef- 
fective relativistic degrees of freedom. On a given scale 
l p , a viscous regime corresponds to R e {l p ) <C 1. Turbu- 
lent flow corresponds to R e (l P ) 3> 1 where the dissipa- 
tive time-scale is much greater than the eddy-turnover 
time-scale r cddy = l p /v = al c /v. We can now define the 
number of eddy-turnover times N oddy given by integrat- 
ing dN cddy = di/r cddy . Since the epoch we are considering 
is in the RD era we have a = T /T ex t x l 2 and 



tl = TT\ — -- 



2f 



(14) 



where we normalised a$ — 1. The radiation temperature 
today is given by the CMB temperature Tq = 2.725 K. 
With a constant v, given in eq. (jlip. we can integrate for 
N cddy from some initial epoch defined by a temperature 
Ti, thus 



N My {T,l c ) 



vmp 




(15) 



The number of eddy-turnover times A cddy > 1 so that the 
magnetic field amplification is not suppressed. 

In the early Universe, before neutrino decoupling, neu- 
trinos have the longest mfp and are thus most efficient 
at transporting momentum and heat. At high tempera- 
tures the shear viscosity due to neutrinos is low and the 



plasma is in a turbulent regime. However, the neutrino 
mfp increases as the Universe expands and cools, leading 
to a viscous regime. Eventually the neutrinos decouple 
at T ~ 2.6 MeV and the fluid flow becomes turbulent 
once again since now the momentum is transported by 
photons which have small shear viscosity at this time. 
Once the temperature decreases below the electron mass 
T < m e — 0.511 MeV, e^ pairs begin to annihilate and 
the photon mfp increases rapidly. The e annihilation 
completes at around T < 20 keV. After this, as the tem- 
perature drops further, photons begin to diffuse followed 
by photon drag and the fluid is in a viscous regime once 
again 
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The diffusion will also damp away the ve- 
locity fluctuations, see for example ref. [3fJ. Therefore, 
before diffusion sets in, there is always a turbulent regime 
if we consider large enough scales. However, we note that 
the eddy-turnover time increases for larger scales. 

Let us consider the epoch of e ± annihilation so that the 
initial temperature is T — m e . At this time, at temper- 
atures ranging m e > T > 20 keV, the comoving photon 
mfp is given by [3l| 



P = a~ l V 

CTT(™pa 



+ n 2 ) 1 / 2 <7Tn pa 



(16) 



where ot = Stto 2 /3m 2 , is the Thomson cross section, 
a s=a 1/137 is the fine structure constant and n palr the 
number density of e pairs given by 31 1 



2m e T 



3/2 



-(-^) ('+!£]• ' 



At lower temperatures 12 — a 1 /aTn e , where the number 
density of free electrons is 



n e = X e n b = X e Q, b I — 

m p \T 



(18) 



and m p is the proton mass. Since the plasma is fully 
ionised at this temperature the ionization fraction is 
X e = 1. 

The SSD amplification is efficient, i.e. not suppressed, 
for scales which have at least one eddy-turnover time 
whilst the kinetic Reynolds number is greater than unity. 
Let us define the final temperature 7/ where turbulence 
is shut off by setting R e (T = T f ) = l |32J|, i.e. 



R e (T = T f ) = 



vIcQ-tTq n pair (T/) 
(37/ 5 5*) T f 



1. 



(19) 



We can rearrange eq. (TT5|l for T = Tf and substitute it 
into eq. (|19l) . from this we can find the largest possible 



scale l c with N oddy (Tf) ~ 1. With g* ~ 3.36, for the tem- 
perature range considered and the typical velocity fluc- 
tuations in eq. (|11|) . we find the scale l c ~ 0.27 pc. This 
scale corresponds to R e (Ti) ~ 10 12 and Tf ~ 19.6 keV. 
This final temperature is close enough to the temperature 



at the end of e completion to justify the approximation 
in eq. (|T6|) . Figure [1] shows the velocity spectrum at the 
final epoch Tf , the scale l c above is clearly well inside the 
horizon where eq. (|11[) applies. 

The magnetic Reynolds number is R m = AiTal c va, 
where the conductivity at T < m e is given by |33| 



a In A 



2T 



3/2 



(20) 



and A = (l/6y/¥a)^yrnf/n e (T/rn e ). We can estimate 
the minimum conductivity in the epoch considered, 
which occurs at T = Tf, and thus the minimum R m . 
With eq. dHJ we find R m ~ 10 15 for l c and T) given 
above, and thus the Prandtl number at this tempera- 
ture is P m = R m /R e ~ 10 15 . This large Prandtl num- 
ber means that we can neglect dissipativc effects due 
to finite conductivity. The large magnetic Reynolds 
number R m ^> i?£f ~ O(10 2 ) [l5j justifies our choice of 
R e (Tf) — 1 as the final state of efficient amplification. 

For the length scale given above, the magnetic field 
growth rate T is much larger than the Hubble rate. 
From eq. (|12[) . at the initial epoch considered, we find 
T/H ~ 10 11 . Due to the exponential growth, the mag- 
netic field will saturate very quickly within only one eddy- 
turnover time. Saturation is given by the equipartition 
between kinetic and magnetic energy. With no net mag- 
netic helicity this occurs when [29( 



(B 2 )^4ire(p + p)(v 2 ) 



16tt 3 £ 
90 



5*TV), 



(21) 



which we evaluate at Tf and e is the saturation effi- 
ciency (see discussion below). We find, with the ki- 
netic energy given by the typical velocity fluctuations 
in eq. ([IT]). B sa 3.7e 1/2 x 10 6 G, whose value today is 
B = B(T /T f ) 2 = 0.35s 1 / 2 nG. Although those strong 
fields are generated at an early epoch in the evolution 
of the Universe, they will not suffer further dynamical 
damping due to the flat velocity spectra and the sur- 
vival of overdamped Alfven waves in the cosmic plasma 



28l l29i . |34j . The generation of a tiny seed field together 
with this efficient amplification via the SSD mechanism 
can easily explain the observed large extragalactic mag- 
netic fields [5[ and have a significant effect on a number 
of astrophysical processes. 



DISCUSSION AND CONCLUSION 

As noted in a number of numerical and analytical 
works, rapid amplification of seed magnetic fields occurs 
due to the turbulent motion of the conducting plasma. 
This small-scale dynamo (SSD) mechanism is believed 
to play a crucial role in the formation of large magnetic 
fields in a number of astronomical settings, from stars to 



paper we have demonstrated that the conditions neces- 
sary for such turbulent amplification arise in the radia- 
tion era before the onset of structure formation. Turbu- 
lent plasma motion arise inevitably from perturbations of 
the gravitational potential; required for structure forma- 
tion and believed to originate from quantum fluctuations 
of suitable fields during a period of cosmic inflation. Seed 
magnetic fields will almost certainly be generated at some 
level in the early Universe through a variety of mech- 
anisms. Such mechanisms include inflation and phase 
transitions. A particularly natural mechanism, proposed 
by Harrison 14 1 and implemented in refs. 10Hl3|. makes 
use of the vorticity generated by density perturbations 
and inevitably leads to the generation of seed fields. The 
generation of seed fields together with the SSD mecha- 
nism for their amplification could explain the observed 
large magnetic fields throughout the Universe. 

The SSD mechanism is very efficient for rotational 
velocity fluctuations, where numerical simulations have 
shown that the saturation efficiency e is close to unity 
24|. The saturation level is lower for irrotational veloc- 



galaxies and the intergalactic medium [lfj, |2l|, |22j . In this 



ity fluctuations e ~ 10~ 3 — 10~ 4 [24j, as we consider here. 
However, these numerical simulations were performed in 
a non-relativistic setting c s «l and at relatively high 
Mach numbers. The required further numerical work to 
investigate the saturation level at low Mach numbers is 
left for future publications. 

Rotational fluid velocity fluctuations v± are generated 
at second order in cosmological perturbations. Indeed, 
the generation of vector modes by density perturbations 
at second order, which contribute to the fluid 3- velocity, 
has been investigated by a number of works [35l438| | . In 
ref. [36(] , the spectrum of vector modes was calculated for 
the radiation era. The power spectrum Vy/A^ grows 
as y on superhorizon scales, peaks when inside the hori- 
zon and then decays as y~ A [36(. When the power spec- 
trum peaks Vv ~ A^. Hence, the maximum rotational 
velocity fluctuation is v± ~ \fPv ~ 10 -9 and could con- 
tribute efficiently to the SSD action. However, the vec- 
tor modes have a x 2_ distribution [36|. It is not clear 
how this non-Gaussianity affects the SSD action, fur- 
ther work in this direction is needed (although authors in 
ref. [331 suggest that breaking the Gaussian assumption 
does not change the qualitative results). Observational 
constraints however allow much greater amount of vec- 
tor mode. In ref. |37| . neglecting gravitational waves, the 
authors showed that V v /V s < 0.012 from WMAP 7-year 
data. If such large vector mode exists, the SSD amplifica- 
tion could be very efficient in the early Universe. Indeed, 
this opens up the intriguing possibility of efficient mag- 
netic field amplification due to vector modes generated 
in non-standard models of inflation. 

Finally, we stress that the main assumptions of the 
Kazantsev model are simplifying assumptions needed to 
obtain an analytical grip on the problem. These assump- 
tions regard the statistical properties of the velocity field. 



The velocity field is taken to be Gaussian, homogeneous 
and isotropic in space and instantaneously correlated in 
time. So far we have not yet calculated the detailed sta- 
tistical properties of the velocity fluctuations. We have 
only seen that the irrotational component has a Gaus- 
sian distribution whereas the rotational component is 
non-Gaussian, and neither are instantaneously correlated 
in time (see for example ref. [40| for related numerical 
work) . For this reason it would be of high interest to de- 
velop an analytical model for more general velocity fields. 
We leave this work for future publications. Nevertheless, 
we expect that below the integral scale of the irrotational 
velocity fluctuation spectrum Kolmogorov type turbu- 
lence is established. But we stress that the Kazantsev 
mechanism will work independently of the type of turbu- 
lence [41|. Hence the early strong magnetisation of the 
Universe seems unavoidable. 
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